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Abstract 

The off-shell two-loop correction to the massive quark propagator in an ar- 
bitrary covariant gauge is calculated and results for the bare and renormalized 
propagator are presented. The calculations were performed by means of a set of 
new generalized recurrence relations proposed recently by one of the authors 
From the position of the pole of the renormalized propagator we obtain the re- 
lationship between the pole mass and the MS mass. This relation confirms the 
known result by Gray et al. 0. The bare amplitudes are given for an arbitrary 
gauge group and for arbitrary space-time dimensions. 
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1 Introduction 



In electroweak precision physics there is ongoing interest to improve the status of 
precise theoretical predictions of the Standard Model. Attempts are being made by 
several groups to provide more complete two-loop calculations. The major obstacles in 
doing so are still lacking methods to perform accurate numerical calculations of two- 
loop integrals. Unlike in the one-loop case integrals cannot be performed analytically 
in general. Recently substantial progress was made in the development of efficient 
algorithms which allow us to calculate two-loop diagrams for arbitrary masses in a 
canonical way [lj . In the present paper we describe a simple application of these tech- 
niques: the calculation of the two-loop quark propagator in Quantum Chromodynamics 
(QCD) with one massive quark and the others massless. Bare expressions are obtained 
for an arbitrary gauge group and for arbitrary space-time dimensions. The knowledge 
of the full off-shell propagator is of particular interest for the discussion of the scheme 
dependence of radiative corrections, since it allows us to relate not only the pole mass 
M to the MS mass m but also to the so called momentum subtraction schemes 0] 
(MOM) and the Euclidean mass definition ||. In lattice QCD a MOM subtraction 
scheme was recently proposed for the non-perturbative operator renormalization with 
matching conditions to perturbation theory far off-shell ||. These considerations are 
of great physical interest in particular in the context of decoupling of heavy particles 
and the effective QCD schemes where it is required to relate QCD in the MS scheme 
for effective theories with different number of (light) flavors Q. 

The two-loop corrections to the electron propagator were first considered in Quan- 
tum Electrodynamics (QED) 0. In QCD so far only the relationship between the pole 
mass and the MS mass was investigated to two-loops 0. A similar on-shell calcula- 
tion, within the dimensional reduction scheme, was given in ||. A detailed discussion 



of the mass renormalization was presented in [H for QED, in |10[ for the electroweak 



Standard Model and in jO]] for grand unified theories. 

In sect. 2, we introduce some notation for the quark propagator. The calculation 
and the results for the renormalized amplitudes are discussed in sect. 3, while sect. 4 is 
devoted to the relationship between the pole mass and the MS mass. Technical details 
and bare amplitudes will be given in a few appendices. 



2 The quark propagator 

The renormalization of the fermion propagator is conveniently discussed by looking at 
the full inverse bare propagator first 

^ 1 (g)=^-m -S(g), (1) 

where S(p) is the one-particle irreducible fermion self-energy. We write the inverse 
bare propagator as 

S F 1 ( ( l) = i B baxe ~ "^bare (2) 

in terms of two dimensionless scalar amplitudes A bare , -Bbare, which are functions of 
the four-momentum square q 2 , the bare mass mo, the bare QCD coupling constant 
a s0 = <7o/(47i") and the bare gauge parameter £ - Off the mass-shell, when q 2 ^ M 2 , 



the bare amplitudes are infrared finite but ultraviolet singular. As usual we start from 
dimensionally regularized bare amplitudes. Renormalization of the fermion propagator 
is performed order by order in perturbation theory by writing the bare mass as a sum 
of the renormalized mass plus a mass counterterm 

m = m + 5m = Z m m, Z m = l+ 6 -^ (3) 

m 

and, correspondingly, for the bare strong coupling constant 

Sex 

a s o = a s + 5a s = Z g a s , Z g = H , (4) 

while for the bare gauge parameter we have 

fo = Z 3 f. (5) 

We define £ such that £ = 1 in the Feynman gauge and Z 3 is the wave-function 
renormalization factor of the gluon field. Dimensional regularization is used with d = 
4 — 2e being the dimension of space-time. The Feynman diagrams are divided by 
27r d / 2 exp(— 'ys)fi~ 2£ per loop with [i the MS renormalization scale, 7 the Euler constant. 
The renormalized fermion propagator is obtained by multiplicative renormalization 
with a suitable wave-function renormalization factor Z 2 

S F ven(q) = -^S F (q) . (6) 

We write the inverse propagator in the form 

rente) = fa B Tcn - 171 A ren . (7) 

Explicitly, we obtain the renormalized amplitudes in the MS scheme by applying the 
following MS renormalization constants [12] : 

+ ("(IT + f + 1t)Ca + tn f + -C F ) - C F , 



„ „ a s /,13 8 \ 1 , . 

' 4ir V 3 3 ') e 

z ™ = 1 * S {3Cf) I + (S) 2 c " ( ( ¥ c * - 2tn > + 1 Cf, 7< 

f 9T 5 , 3 , 1 
- (— C A - -tn f + -C F )- 
v 12 3 i A e 

For SU(N C ) we have C F = (N 2 C - 1)/(2JV C ), C A = iV c , t = 1/2. Since the latter 
coefficient does not depend on the number of colors N c its value will be inserted in 
most of the formulae presented below. 



3 Renormalized Amplitudes 



There are six diagrams contributing to S(g) at the two-loop level (see Fig. 1). 





(c) 






Fig. 1. Fermion self-energy diagrams to two- loops. 
Solid lines denote quarks, wavy lines gluons and 
dashed lines Fadeev-Popov ghosts. 



In our calculation we consider the linear covariant gauge with arbitrary gauge pa- 
rameter £. UV and IR singularities are dealt with by dimensional regularization and 
cancel in the observable quantities. A major problem in performing higher loop calcu- 
lations with arbitrary masses is to find an appropriate basis of integrals which allows 
us to present results in a compact form. In most cases integrals can be reduced to 
a set of master integrals by using partial fraction decomposition, by differentiation 
and by integration by parts. The remaining problem of integrals with irreducible nu- 
merators was solved recently in [|l3j. Such integrals may be expressed in terms of 



scalar integrals with dimensions shifted by multiples of 2 which then are reduced again 
to integrals of dimension d = 4 — 2e (generic). The resulting complete set of rela- 
tions between integrals, a set of generalized recurrence relations, was considered in [|I|] 
and has been implemented in a computer program written in FORM ||14|| . All dia- 



grams finally could be expressed in terms of 5 two-loop integrals J3, J m (m 2 , m 2 , m 2 ), 
Jii2(fn 2 , fn 2 , m 2 ), J m (0, 0, m 2 ), Jn 2 (0,0,m 2 ) and 3 products of the one-loop integral 
structures Gn(0,m 2 ) and Goi(0,m 2 ), where 



h{q 2 /m 2 ) = -q- 
j 1 2 2 2 



d d ki d d k 2 



ii d k\{k\ — m 2 )({ki — q) 2 — m 2 )(A; 2 — q) 2 ({ki — k 2 ) 2 — m 2 ) ' 

d d h d d k 2 



n d {k 2 - m 2 ) a ((k 2 - q) 2 - mlf{{ki - k 2 ) 2 - m 2 )^ ' 
d d k x 1 
7T d / 2 {k 2 - m?) Q ((fci - q) 2 - m 2 y ' 



(9) 



The two-loop integrals Jm, Jn 2 and the one- loop integral Gu with zero masses can be 
expressed in terms of hypergeometric functions |l5j as follows 



J m (0,0,m 2 ) 
J 112 (0,0,m 2 ) 
Gn(0,m 2 ) = 



-i 2 (m 2 ^ 



r(i-^)r(^-i)r(4-d) a F 1 



d Q 
2' ° 
d . 
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— d ; 


m 2 


^ 4 
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^m 2 ) d ' 2 - 2 Y{l~ d -) 2 F l 



1,2 



2 ' 
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(10) 



In appendix A their expansions up to order 0(e 2 ) are given. Goi is elementary: 

Goi(0,m 2 ) = -^(mY/ 2 - 1 r(l-^). 

Finally, Jm and J\\ 2 with equal masses may be written in the form 

3 1 fz 9 



J m (m 2 , m 2 , m 2 ) = m 2 " 4 T 2 (l + s] 
J U2 (m 2 , m 2 , m 2 ) = m" 4£ r 2 (l + e) 



2e 2 e V4 



Js 2 



— o V 3 h - 

2e 2 2e 6 3 V ; 



1 7/ 



(11) 



where 2 = q 2 /m 2 and J^(2) = zd/dzJ^z). The numerical evaluation of I3, J3 and J3 
will be discussed in appendix B. Let us mention that results from individual diagrams 
are rather lengthy such that we must refrain from presenting individual contributions. 

The amplitudes in (0) are calculated from the diagrams of Fig. 1 where diagram 
(b) was taken with n/ — 1 massless fermions and one massive one in the loop. We find 



A ren = l + |^ F {4 + 2£+(£ + 3)(^ln2 + J M )j 



Air) 



1 260 103 4 
1 1 h-z + 8f + 



1 z e-{i-^+e 



U 2 {z) 



1 - 31z + I82 2 - £z(9 - 5z) - £ 2 (1 - z 2 )) In 2 + 2(f + 2)(f + 3) J M 
2(l+£) 8(2 + 3) 8(2-9) 

- 4 (l - 52 + 82 2 + 6^ 2 - £ 2 (1 + 2)) In 2 2 
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+ C f Ca 
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J 3 + 



32 
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122 22^ 

+ 5 89 + 1 ^ + 3e2 L M + 1(^ 2 + 3£ + 22)(2ln 2 + |l„)L, 



+ C F t 
52 
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- 2n f Ll + ^(314 - 95z + 15z 2 + 2z 3 ) - 
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(4z + 6z£-z£ 2 ) (C 2 + Li 2 (z)) 



+ ^- (8z 2 + 183,2 - 1140 + 48£(8 + £)) 
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_ _L (7z + 4)n/ _ 2n/ L, - 2^((1 + z)n, + ^(z 2 - 17 z - 14)) lnz] } , 



(13) 



where z 



1 — z and L, 



defined by the integrals 



ln(/i 2 /m 2 ). The polylogarithms Li„(z) and S ntP (z) are 



Sn,p(z) 



_ 1 )n+p-l ,1 In™- 1 ^) l n P(l _ a^) 



(n — 1)! p! Vo 



X 



cfe ; Li n (z) = S n _ ltl (z) 



and Cn = C( n ) — J2i^i^/n l is the Riemann zeta function, with values ( 2 = n 2 /Q, 
C 3 = 1.202057..., C 4 = vr 4 /90 etc.. 

We observe that in order to evaluate the on-shell values for A rcn and B rcn (z — > 0) 
we need to know the on-shell value of Iz{z) and the expansions of Jz{z) up to order 



z 2 and of J^{z) up to order z. The on-shell value of Iz{z) (at z = 1) is known 
and is given for completeness in appendix B. The expansions of the integrals J m and 
J112, and the corresponding expansions of J%(z) and J^(z), can be obtained by the 
second-order differential equation 

2(q 2 + m 2 )(q 2 + 9m 2 )m 2 j!^ - ((d - 4)g 4 + 10(3d - 10)q 2 m 2 + 9(5d - 16)m 4 ) ^ 

48q 2 m 2d-& 



+ 3(3d -8)(d- 3)(q 2 + 3m 2 ) J 3 



(d - Af 



where primes denote differentiation w.r.t. m 2 , which was derived in [17|. As a result 
we find 

, x 59 3 .5 3.0 

/,w = -f - <t - §« • 

In appendices D and E we also present the small and large momentum expansions 
of the renormalized amplitudes A ren and B ven . For details concerning the expansions 
of h(z) and J${z) we refer to appendix B. 



4 Connection between pole mass and MS mass 

Due to the confinement of QCD, quark masses in general do not have an unambiguous 
simple physical meaning. This is especially true for the masses of the light quarks 
which are usually parametrized by MS masses, the MS renormalized Lagrangian mass 
parameters. Physical values may be attached to them as current quark masses which 
cause the observed chiral symmetry breaking (see e.g. For a very heavy quark, 



like the top quark, the situation is different. The top quark is so unstable that it decays 
by weak interactions before strong interactions come into play and form bound states. 
In this case the complex pole mass gains a physical meaning as it manifests itself in a 
resonance peak in physical cross sections. In the following we calculate the pole mass 
at the two-loop level in perturbative QCD. We should mention here that the true pole 
mass of full QCD is expected to differ by non-perturbative effects from the pole mass 
obtained in perturbative QCD. 

The relationship between the pole mass and the MS mass at two-loop order was 
first obtained in 0. The two- loop correction turns out to be rather large. This may 
have physical reasons jl9H , but there could also be a problem with the procedure of 



calculation adopted in [[J. In [0 prior to the e-expansion and the renormalization, the 
external momentum was taken to be on-shell and recurrence relations were applied di- 
rectly to the on-shell integrals. Note that the corresponding interchange of integrations 
with taking the on-shell limit could cause problems due to possible infrared singular- 
ities. One expects that such an on-shell algorithm can be applied for infrared-safe 
quantities. However, since we are lacking a rigorous proof that the recurrence relations 
for the on-shell integrals yield the correct answer, it is highly desirable to find the 
connection between the pole mass and the MS mass by means of a different method 



which avoids possible ambiguities. This is possible by determining the position of the 
pole of the MS renormalized propagator. Thus, setting $ = M in and taking the 
limit q 2 — > M 2 in A ren and B ren , we obtain the transcendental equation 

lim [mA ren (g 2 , m 2 ) - MB ren {q 2 , m 2 )} = (14) 

q 2 ^M 2 

for the pole mass M. Its iterative solution yields the desired relation 



M = ™( 1+C 'IS +C HS/ + -| (15) 

with 

d = C F (4 + 3L), (16) 
c 2 = - 8C 2 - 4J 3 (1) + + ^L 2 



'71 26 
C l An f [- + 8( 2 + -L + 2L 



/ 121 27 9 \ 

+ C 2 + 30C 2 + 8/ 3 (l) + -L + -L 2 j - 12CVt (1 - 2( 2 ) . (17) 

where ^(1) is given in appendix B and L = ln(/i 2 /M 2 ). Setting /i 2 = M 2 we arrive at 
the relation already known from [fj]. 

Note that as long as we consider QCD corrections only, the location of the pole 
of the quark propagator M is real. Only after switching on the weak interaction the 
quarks become unstable and the pole moves off the real axis. 

Our result confirms the validity of the on-shell approach, i.e. the possibility of 
setting the external momentum squared on-shell and applying the on-shell recurrence 
relations from the very beginning. This check thus enhances our confidence that it is 
possible to derive the same relation on the three-loop level by starting from the on-shell 
recurrence relations, which for this task evidently is much simpler than its "indirect" 
evaluation from the renormalized propagator. 



After completion of this paper we learned about related work [2(J where the 2-loop 
relation between the quark mass in the MS and MOM schemes has been obtained in 
the NNLO approximation. The result |20j corresponds to the limit q 2 /m 2 — > oo of the 



complete 2-loop massive quark propagator which has been presented here. 
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Appendix A 

Here the two-loop integrals Jm(0, 0, m 2 ), Jn 2 (0, 0, m 2 ) and G\\ of fllQ|) are expanded 
in e. This is the form which is needed to obtain the renormalized amplitudes A and B. 

e 2 ^m 4 - 2 J m (0, 0, m 2 ) = -± + ^ - 3 + H . - ^2 + In z" - Li 2 (z) 

/ 15 115 3(z-6), 4 A l-z 2 /13 , \, 

+ e 1 z + — -t 2 + -C3 H Inz Inz 

^ 4 16 4 ^ 2T 3^ T z ^4 y 

2z 2 -6z-l T ... T . , , ._ , ,\ 2 /21 865 , 39,. 
+ Yz Li 3 (z)- 4^,2^)1 + e 2 {- + — z -% 2 + —zC,2 

+ (4 - -z)C 3 - — C4 + [ — + -( 2 - — Inz + - In 2 z) Inz 



V 8 

■ 3Lil(z) - Li 4 (z) - ^(z 2 + 6z - 2)5 1)2 (z) + ^(2z 2 - 6z - 1) Li 3 (z) 

" (6 + § " y z + 3C2) Li2(z) + 3(1 ~ z z2) ln ^ Li2(z) _ 165i < 3(z) + 852 > 2(z) 



(18) 



where S np (z) is the generalized Nielsen polylogarithm |21 



e 2 ^m 4e - 4 J 112 (0,0,m 2 ) = __L - 1 + I - ^2 + f Inz - Li 2 (z) 
+ e hr ~ 9^ 2 + — ln z In 2 ^ + -Cs - - Li 2 (z) - Li 3 (z) 

-4S lj2 (z)) + e 2 f f + ^2 - ?C4 + - (19 + 3C 2 ) In z - 10- In 2 z 
V 2 2 8 z z 

8z 4 1 5 z 

+ — In 3 z + -Cs - - Li 3 (z) + Li 2 (z)(6 3C 2 + 6- Inz) 

6z 6 z z z 

+{- z - 12)S 1)2 (z) -3Li 2 (z) - Li 4 (z) - 165i, 3 (z) + 85 2)2 (z)) . 



(19) 



e 2 7£m 4 £ -2 Gii ^ m 2 )G()i ^ m 2 



1 1 

'T 2 ~ ~£ 



3+ - Inz ) - 7-C2 + -(-3 Inz + In 2 z + Li 2 (; 

z I z V 



+ £ 



-15-3C 2 + ^C3 + -(-7-C2 + 3 lnz- ^ln 2 z)lnz 
3 z V 3 



+ - 1 (3 - 2 ln zj Li 2 (z) + Li 3 (z) - 2Si )2 (z) 



2 Z 

+ £ - 

z 



J (-31 - 7C 2 + 2C 3 - 7 -U) ~ (l5 + 3C 2 - ) In : 



+ (7 + C 2 )ln 2 z- 21n 3 z + - In 4 z + (7 + C2 - 61nz + 21n 2 z) Li 2 (z) 
+ (3-21nz)(Li 3 (z)-2S' 1]2 (z)) + Li 4 (z) + 4S lj3 (z) - 25 2 , 2 (z) 



(20) 



The last formula can as well be used for the e-expansion of the one-loop scalar propa- 
gator integral with arbitrary masses and momentum (see, for example, p2|). 



A surprising fact is that in the final renormalized result of the quark propagator 
given in sect. 3 various structures like £ 3 , £ 4 , Lil(x), Li 3 (x), Li 4 (x) and all generalized 
Nielsen polylogarithms, except Si^, have canceled. 

Appendix B 

In this appendix we present details about the calculation of the functions I 3 (z) and 
J 3 (z) entering the self-energy function given in sect. 3. We can give one-dimensional 
integral representations which are, however, not the most convenient forms for a nu- 
merical evaluation as will be seen. In particular it turns out that the expansions w.r.t. 
small and large q 2 are very efficient even in the time-like region (on the cut) of the 
functions. For this reason we give a full account of these expansions, mainly due 



to p|. 



An integral representation for J 3 is given by the following expression [|16| |17 

2 2 2 2 

w w — m w — m 

„ — In 5 — 

m 



q 2 I 3 (q 2 /m 2 



00 4w dw 



m w 2 — m 2 



In 



w 2 + q 2 



m" 



+ 



4 dw 



31n — 

m 



w 3w 2 -3m 2 + q 2 , W+ + W. 



w 



W+W- 



/2m (w 2 — Am 2 ) 1 / 2 
with W± = ((w± m) 2 + q 2 ) 1 / 2 . 

We recall that for q 2 = m 2 (z = q 2 /m 2 = 1) the value of I 3 is WE 



In ' ^ ' — 1 In , (21) 

W+-W-J 2m' v ' 



-C 3 -6C 2 ln2. 



(22) 



The small q 2 expansion is obtained from |T6| in the following manner: with the 
Ansatz 

Js(z) = - £ z n ^- , (23) 

n=l 



n 



d 3 (n) being expansion coefficients and rf 3 (l) = — C2 + y^ 2 w ^ n ^2 derived from the 
maximum value of Clausen's integral, 

S-2 



JLci(-) = 0.2604341376 . . 
9y 3 3 



(24) 



The higher order terms are obtained by solving the following differential equa- 
tion JL6| (the notation is adopted from [[HJ) 



where 



D 3 D 2 I 3 (z) = D 3 D 2 I a (z) + (q 2 - 9m 2 ) 2 - 2<K 2 m 2 (q 2 + 3m 



d 



(25) 



D 2 = (q 2 -m 2 )-, 
D 3 = 6m 2 (q 2 + 3m 2 ) + (q 2 - 9m 2 ) 



d 



(q — m ){q —9m)-— — 8m 
dq z 



2 d 
^ dq 2 



and 




n-l 

(26) 



Inserting fl23|) into (p5|) , equating equal powers of z and solving iteratively for the 
coefficients d%{n) is very conveniently done, e.g., with the help of a REDUCE |23| pro- 
gram, which can also perform the final numerical evaluation. To calculate J3 on its cut, 



the method of mapping and Pade approximation is applied as introduced in j24j . Since 
any number of Taylor coefficients is easily obtained for J 3 , it can be calculated with 
practically any required precision even on the cut. Thus as an example for z=100 with 
100 Taylor coefficients we obtain J 3 = 7.77982600 — 0.779167952, where the precision is 
estimated from the convergence of the Pade approximants. For comparison: the inte- 
gral representation (|2T| ) is not even applicable on the cut and in the space-like region 
for such large q 2 the two contributions mainly cancel, which seems not the proper way 
to evaluate ^3 numerically. 

As described in [16|] as well, for large q 2 the large momentum expansion can be 
used: with five terms 

(h)^) = 6C3 + (2L 2 Z + 6L Z + 6)/z + (4L 2 Z + \h z - ^)/z 2 + 

,29 r2 46 , 257 w , , r2 1957 r 3613 w 4 

( — L 2 L z z 3 + (32L 2 L z ) z A + 

v 3 z 3 18 ;/ v 2 24 96 ;/ 

,677 . 30907 T 103577, . «, 

( K L z )/z 5 , 

V 5 2 75 1000 ;/ ' 

and L z — In z — m we obtain the same numerical result as above except for the last 
decimal in the imaginary part (i.e. Im (73) (100) = —0.77916794). 

We also mention that from the small and large momentum expansions of Iz{z) a 



3-loop integral D 3 was calculated [25 



The situation is similar for J 3 . We found the following integral representation 

n 13 (3 + ^)(l-^) , 

J 3 (z) = -9 + -z+ K - g ln(l - z) 

, n s f 1 dy[-3 + y(3-z)-3y 2 ] ( 1 - y + y 2 - zy \ 

+{1 - z) i ^yJ^W 1 1 7K )^ (27) 

with A = y A — 2(1 — z)y 3 + [z 2 — 6z + 3)y 2 — 2(1 — z)y + 1. However, this representation 
is only applicable for z < 1 . Recall that the first threshold starts at z=9 ! Therefore, 
the small and large momentum expansions will be much more useful again. For J${z) 



we have the special values p6 | 



Writing 



j 1 59 t / \ 45 8tt 
A(l) = - T , *(9) = T -^ 



n=0 



the Taylor coefficients j 3 (n) of J 3 are obtained recursively in the following manner 



21 27 S 111 1 Q 1 

J 3 (0) = -- + y 5 2 , j 3 (l) = g - 35 2 , j 3 (2) = — - -S 2 , j 3 (3) = — - -S 2 . (29) 



Introducing c 3 (j) = —j\(J + l)!j 3 (j), j > 2, we have the recursion JT7 

c 3 



(J + 1) = I ((10j 2 - 4)c 3 (j) - (j - 2)(j - l)j(j + l)c 3 (j - 1)) . (30) 

Thus the coefficients are much easier to calculate than in the case of J 3 . Furthermore, 
the precision of the Pade approximants is much higher, i.e. already with 33 coefficients 
one obtains for z = 100 a precision of 10 decimals. The reason for this is of course that 
the threshold is much higher in this case. 



The large q 2 expansion of J 3 is given in [jlT | . The first six terms (again up to z 5 ) 
give the contribution 

(^) asy (*) = - f)z + \l? z - H _ {3 Ll + \ Lz - \)/z 

- m - QL Z - |)/z 2 - (9Ll - 23L Z - ^)/z 3 - (36L 2 - 107L, - ^)/z 4 

-(171L-^-^. (31) 

For 2=100 this yields also a precision of 10 decimals, i.e. the same precision as the 
small q 2 expansion with 33 coefficients. 

What concerns J 3 , the Taylor series as well as the large momentum expansion are 
easily obtained by differentiation from the corresponding representations of J 3 . The 
precision with the same number of coefficients as above is here 9 decimals for the small 
q 2 expansion and the asymptotic expansion. 

Finally we mention that also the Taylor coefficients can be calculated analytically 
and stored which allows us to write a FORTRAN program in terms of the multiple 
precision version of D.H. Bailey ^7|. Then, calculating these functions from the small 



momentum expansion, for J 3 , J 3 and J 3 only the Pade's have to be performed, which 
is very little time consuming. For the numerical evaluation of Si t 2 we refer the reader 



to |2T . 



Appendix C 

In the following we present the bare amplitudes A hare and B hare in terms of the 
basic integrals of sect. 3. For these we introduce the following notation: 

7i = Gii(0,m 2 )Gn(0,m 2 ), 
h = -4gi,i(0, m 2 )G Q1 (0,m 2 ), 
I 3 see (|), 

T T I 2 2 2\ 

vnr 

h = Jii2(m 2 ,m 2 ,m 2 ), 



h = -^Jm(0, 0,m 2 ), 



m 2 

h = Jn 2 (0,0,m 2 ), 

/8 = ^rG'oi(0,m 2 )G'oi(0,m 2 ) 

The general structure of the bare amplitudes reads 

C F (C A ai C F a 2 \ T C F ( ( „ 1„ N , ta A (3d-8) 
A = — Id— 4 + —) h + —z [ {Cf ~ 2° A)a3 + ~W^W 
(d-2)C F [C A a 5 (l-z)(d-3)t C F a 7 \ T 
+ (l-z)(d-3) [a~l + d^2 a6 + —) 12 

+ 9£(c F - 1 -C A )a 8 I 3 + ^-((C F - 1 -C A )a 9 + ta ™ V > 



2 



z\ r 2V°° 1-zV 2 " 3(l-z) 

(3d-8)C F ( C A a u ta 12 \ 
+ (d-4)(l-s) W = 4) + 3d~8 + CFai3 J h 
Cf I C A a u Cfciiq\ 
+ ( d -4)(l-*) {(d-4)(d-3) + tal5 + 1T3 ) h 

(d - 2) 2 C F ( C A a 17 iois CV^ig \ . /o 9 \ 

(l-z)(d-3) \d-A + 3(l-z) 2 (d-5) ^ J 8 1 J 



R _ C F /C7 A 6i n ,\ T (d-2)C F (C A b 3 t(l-z)(d-3)b 4 
ZB ~— [d^A + Cph2 ) h + (d-3)(l-z) V ~T~4 + 15(d-2) 

+ ^ ~ 2^ + 15(1 - ^) 2 J /5 + —z \2{d - 4) 2 (d - 6) + d^6 + d^4 

/ Ca^u ^15 C F biQ \ 

T~^\(d-3)(d-A) 2 (d-Q) d^6 (d-4)(d-3)/ 7 

+ (d - 3)(1 - z) {in + 15(d - 5)(1 - ,) 2 + CFbl9(d " 2) J /8 ' 

where the coefficients Oj and bi ( i=l,2, . . . ) are polynomials in the space-time dimen- 
sion d, the gauge parameter £ and z = q 2 /m 2 . 



ai = (d - 4)(d 2 - 4d + 2) - (d - 2)(d 2 - 18d + 46)z - 2(d - 3)(d - 4)£ - 2(d - 3)(d - 6)£z 

a 2 = (d - 2) 2 + (-4 - 4d 2 + Yld)z + (12 - d 2 )z 2 + 4(d - l)(d - 4)£z + 4(d - l)(d - 2)£z 2 

-(d - 2)(d - 6)£V + (-2d 2 - 20 + 12d)£ 2 z - (d - 2) 2 £ 2 

a 3 = 8d(3d - 8) 

a 4 = -48d + 264 + (-96d + 272)z + (16d-24)z 2 

a 5 = -d 3 - 86 + 47d-3d 2 + (2 + d 3 -9d 2 + 17d)z + (3d-10)(d-5)£ + (d-2)(d-5)£z 

a 6 = 8/3(d- l)(d- 2) 



a 7 = _2(d-2)(d-3)-2d(-13d + 22 + 2d 2 )z + (-2d + 12)z 2 

-2(d - l)(5d - 16)£z - 2(d - l)(d - 2)£z 2 + 2(d - 3)(d - 4)£ 2 z + 2(d - 2)(d - 3)£ 2 
o 8 = 24 + 2d 2 - 12d - (2d 2 + 8 - \2d)z 
ag = 8zd — 72d 

do = -2160 + 432d + (-2352 + 816d)z + (432 - 240d)z 2 + (-16 + 16d)z 3 

o n = - (d - 4) (d 2 - 9d + 16)£ 2 + (2d 3 - 8 - 18d 2 + 40d)£ - d 3 - 7d 2 + 60d - 104 

a 12 = -8(d-2)(3d-8)(n/- 1) 

o 13 = 2(d-3)(d-6)£ 2 -8(d-l)(d-4)£-4-14d + 6d 2 

aw = -216d 2 + 628d + 26d 3 - 640 + 2(d - 2) (d - 4) (d 2 - 3d + l)z + (332d - 4d 4 + 50d 3 

-128 - 210d 2 )£ - 2(d - 5)(d - 2)(d - 4)£z + 2(d - 4)(d - 3)(d 2 - 9d + 16)£ 2 
oi5 = 32(n/-l)(d-2) 

oie = -144 + 68d 2 -16d-16d 3 -4(d-2)(d-4)(2d-3)z 

+4(7d - 22)(d - l)(d - 4)£ + 4(d - l)(d - 2)(d - 4)£z - 8(d - 6)(d - 3) 2 £ 2 
an = (-d + 5)£ + (d+l)(2d-7) 

o 18 = -8d 2 + 204 - 4d+(80d 2 -616d + 1112)z-4(d-l)(2d-9)z 2 
a 19 = d - 3 + (2d 2 - lOd + A)z + (4d - 4)£z + (3 - d)f 



b 1 = -2(d-2)(d-3)^ 2 -2(d-3)(d-4)^ + (d-2)(d 2 -10d + 22)z 2 

+ (56 - 34d + 8d 2 -d 3 )z 
b 2 = -3(d-2) 2 + (36-12d)z-(d-2)(d-6)z 2 + 4(d-l)(d-2)£ + 4(d-l)(d-4)£z 

-(d - 2) 2 fz 2 + (-2d 2 - 20 + 12d)fz - (d - 2)(d - 6)£ 2 

6 3 = _56 + 24d-2d 2 + (d 3 + 6-5d 2 + 5d)z-(d-2)(d 2 -9d+19)z 2 

+ (d - 2)(d - 4)£ + (d - 3)(3d - 10)£z + (d - 2)£z 2 

64 = 2(d-2)(d-4) -4(d-2)(d-4)z + 2(d-2)(d-4)z 2 

65 = -12-10d + 4d 2 + (-20 + 16d-2d 2 )z + 2(d-2)(d-5)z 2 

+ (d - l)(2d 2 - 15d + 26)£ + 2(d - l)(d - 3)(d - 4)£z + (d - l)(d - 2)£z 2 
+2(d - 2)(d - 3)fz + 2(d - 3)(d - 4)£ 2 

6 6 = -2d(d-6) + 2(d-2)(d-4)z 

b 7 = -2(7d-18)(d-6)-10(d-2) 2 z 

b 8 = 36(7d-20)(d- 3) + (-4680d + 576d 2 + 8336)z + (2288 + 408d 2 - 1928d)z 2 

+ (472d - 592 - 96d 2 )z 3 + 4(3d - 8)(d - 3)z 4 
69 = -216 + 36d+(32d-48)z + (-4d + 8)z 2 

b w = -756d + 2160 + (8832 - 1644d)z + (4d - 16)/ + (2160 - 1032d)z 2 

+ (-1072 + 424d)z 3 + (224 - 68d)z 4 
6 n = -3168d-8d 5 + 2624d 2 + 132d 4 - 874d 3 + 768 -(d-2)(d 4 + 5d 3 -92d 2 + 312d 

-320)z + 4(3d - 8)(d - 4)(d 3 - 9d 2 + 25d - 18)£ 

+2(3d - 8)(d - 2)(d 3 - 7d 2 + 4d + 28)£z - (3d - 8)(d - 2)(d - 4)(d 2 - 9d + 16) fz 
-2(3d - 8)(d - 3)(d - 4)(d 2 - 9d + 16)£ 2 
b 12 = 16(d-2)(d-3)(n/-l) + 8(d-2) 2 (n/-l)2 



6 13 = 2(d-4)(5d 2 -16d + 6)-2(2d-7)(d-2)(d-4)z-2(d-l)(3d-10)(3d-8)£ 

-2(d - l)(d - 2)(3d - 8)£z + 2(3d - 8)(d - 2)(d - 3)£ 2 z + 4(3d - 8)(d - 3) 2 £ 2 

614 = 2160d - 86d 4 - 1792d 2 - 384 + 5d 5 + 585d 3 + (608d 2 + 960 - 109d 3 + d 5 + 2d 4 

-1304d)z - 2(d - 2)(d - 4)(d - 6) (2d - 7)z 2 - 2(3d - 8)(d - 4)(d 3 - 9d 2 + 25d - 18)f 

-2(3d - 8) (d - 3) (d 3 - 8d 2 + 14d + 4)£z - 2(d - 6) (d - 2) (d - 4)£z 2 

+ (3d - 8)(d - 3)(d - 4)(d 2 - 9d + 16) fz + (3d - 8)(d - 3)(d - 4)(d 2 - 9d + 16)£ 2 

6 15 = -16(n/ - l)(d- 2) - 16(n/ - l)(d- 2)z 

6 16 = _4(d-4)(3d 2 -lld + 3)+4(d-4)(d 2 -9d+13)z + 8(d-2)(d-4)z 2 

+2(d - l)(3d - 10)(3d - 8)£ + 16(d - l)(d - 3) 2 £z - 2(d - l)(d - 2)(d - 4)£z 2 

-4(3d - 8)(d - 3) 2 fz - 4(3d - 8)(d - 3) 2 £ 2 
b 17 = _(d_2)(d 2 -10d + 28)-(d-2)(d 2 -5d+5)z-(d-2)(d-4)C + (-d + 2)^ 
bis = -2(d-5)(d-3)(d-4) + (-72d 3 - 4044d + 988d 2 + 4664)z 

+ (-3676d + 3392 - 140d 3 + 1288d 2 )z 2 + (-212d 2 + 24d 3 + 612d - 616)z 3 

-2(d-5)(d-3)(d-4)z 4 
b 19 = _i5 + d+(2d-2)z-(d-l)(2d-7)^ + (-d+l)^+(3-d)^ 2 



Appendix D 



From the small q 2 expansion, as described in appendix B, we get the following 
expressions in the limit z = q 2 /m 2 — > for the renormalized amplitudes A rcn and B ren . 



A = C F ^-{l + 3L, + f (1 + LfA + z{\ + \0 + z 2 ( l -+ 

- ■ • 5 81 _ r , 313 7 . 3 .9. T o , 11 3 . 1 . v 
■24 + IT* + " + M ~ + ^ + ^ } + L ' ( ~ + ^ + ^ } 



47T 

+(£) 2 {c,C F 



3. l.o, /403 5 . 5.9 9„ , ,11 3. 1 9, ... 1 



24 



2^ 4 

1 

M" 2 r f 



4" 4 
1 

2^ 



+c 2 (i - ^ + ^e 2 ) + + -a + ^ 2 - 7^2 + w{- + 7 e + -e) + c 2 (i + ^ 

n/313 2 A 1 9 1 „ r ,11 1, L 2 > 1, A 
+ ^ ("36 + 3^ + 3^ " 4 * 2 + + 4^ + 12 & + 6^ 2 



+ C 2 F [-19 - 10$ + e 2 + 81S 2 + L M (-15 - 2£ + £ 2 ) + L 2 (- + 3£ + -£ 2 ) + C 2 (l + 6£ - £ 2 ) 

/ 57 , 1,9 9 ,91,, \ 9/ 295 1 1 9 1 

+K"T + ^ + 4^ + 2 S * + L ^-2 + 2* > + 3C2 ) + ^ (-36" + 2^ + 12* + 2* 2 

20 



-32 



20 



+4C 2 (1 - n/) + z(-60 + 252^2 - y«/ - 2^^) 
+, 2 (-^ + 352S 2 - fn, - Hn^)] } + 0(z 3 ) 



n/ + 1625 2 - —rifLf, - 2n/L^ 



(34) 



3 4 



3. . I.,, . /71 , 13, , 1„ 2 „_ , CT , 1 . 2r . ,11 1 



+C 2 (3 - -i + -e) + z{- + -i + -£ 2 - 155 2 + £L„ + -£% + c 2 (- + -e 



+z 
+C 2 F 

+z 



,/1409 5 A 17, 

' 1 — f h — e 

V 432 6^ 48 ^ 

1 \e+^s 2+ L,{-\-^ + \e 



1 QQ 3 1 1 

^ 2 + M-^ + -e 2 ) + c 2 (i + ge) 



^ 2 + 305 2 + L M (-4e + ^ 2 ) 



o " 8C 
49 11 

12 ~ "6 

+^(-3e + ^ 2 )-C2) 

2 / 1628 1 700 x 11 ^, 3 , 



c 2 + * 



599 9 



133 



216 



24 - -n/ + IO8S2 - 2n/L^ + z( — — - -n f + I68S2) 



(35) 



Appendix E 



Using the asymptotic formulae discussed in appendix B, we obtain the following 
expressions in the limit z = q 2 /m 2 — > 00 for the renormalized amplitudes A rcn and B ien 
(L q = ln(-g 2 /^)) 



A30 = -7- Cf 

47T 



(4 - 3L, + {(2 - L,)) + i(l + L, + L„)(3 + - (3 + f ) 



1531 



24 



- 3(7 + 0C3 + 10f + ^ 2 - + + ^ 2 )L 9 



11 



1 



1 /349 



181 



z V 12 

,„ 3 A 1 2 . r2 .313 7 
■(11 + ^+2^2 + (l2- + 2 



+ + 7* + 3 ( 3 + OCa + (^2" + ^ + 



e+^ 2 -(H + ^+^ 2 )^) 



+ 



14 1 
3 2 



+ -( " + (23 + 3? + ^ 2 )L, + J(9 + ()L 2 q 



3 
4 1 



,35 9. 1 . 2 3 3 . 



+ t 



+( y + ^ + ^ + + 0^)^ + 4(9 + 

1 / 8 20 20 

+ ~ z \ 3 LgUf ~ 2ALq + 4L2qUf ~ Y n/ " (24 " 4L,?n/ + Y nf)L ^ 

+ J_ (^48 - 4(n/ + 9)L g + Jn/ - 2(18 + n/)L^ 



, 52 3 2 r2 , 



+ C F 



13 + 12C 3 + 8£ + e - 2(i 2 + 5£ + 6)L q + (H + 3 £ + ^ 2 )L 2 



+ i (-16 + 4£ + 4£ 2 + 12C 3 + 6(3 + £)L q - (£ + 3) 2 L 2 
+(27 + 12$ + f- t 2 L q + 9L,)L„ + 6(3 + 0^) + ^ (~f - y£ + ^ 2 
_( 3 _ 6 £ - 4^ 2 )L, - 6L 2 - + 3£ - ^ 2 + 12L 9 )L M - 6L 2 )] } + o(l). 



(36) 



z z 2V 2 
41 13 9 

T + l2 



25 7 3 

f ^ + ^ 2 - 3(1 + OCa " ( T + 2^ + t/) L « + 7(3 + 



1 

4 l 



1/31 Q 

+ ~z (T + ^ + e 2 + 3«-3)C3-f(3 + OA 



1 / 27 19 3 ~ _ ,17 7, 1 », r ,9 3 A 1 2n . 



+ t 



{2L q - 7 -)n } + - (18 - 12L 9 - 4n f - 12L„) 



1 /332 8 8 \ 



+ c F 



L q e + \L q + X -eL\ + I (-12 + 10* + 4£ 2 - 2£(6 + OA 



+ 4 f- + ^* + ^ - 24Cs + (19* + 3e - l)L q - *(9 + 0^2 



19* - -* 2 + e(6 + t)L q )L^ - 3£L 2 



+ °(72)- 



(37) 
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